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Abstract 

Food webs are networks describing who is eating whom in an ecological community. By 
now it is clear that many aspects of food-web structure are reproducible across diverse 
habitats, yet little is known about the driving force behind this structure. Evolutionary and 
population dynamical mechanisms have been considered. We propose a model for the 
evolutionary dynamics of food- web topology and show that it accurately reproduces observed 
food-web characteristic in the steady state. It is based on the observation that most 
consumers are larger than their resource species and the hypothesis that speciation and 
extinction rates decrease with increasing body mass. Results give strong support to the 
evolutionary hypothesis. 
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1 Introduction 



The complex networks of feeding relations in ecolog ical communities, food webs, have fascinated 
researchers for over t hirty years llCoh en et al lll990l). U sing carefully collected empirical data, 
previous indications l|Cohen et alT 199fl iHavenst 1 199^) could be confirmed that several aspects of 



islands, deserts, and lakes dCamacho et al. 


2002bHMilo et all l2002HGarlaschelli et all 12003: 


IWilliams and Martined 1200(4 ICattin et al 


2004). Most striking, perhaps, were findings that 



certain descriptive model food webs, constructed by following just a few simple rules for 
connecting the elements of an abstra ct species pool, can quantitatively reproduce large sets of 
basic statistics of f ood-web topology l|Cohen et all Il990l IWilliams and Martined . l200fi 
ICattin et all [2004) . These models do not yet expose the dynamical processes or mechanisms by 
which the structures they describe are formed. But they set a standard by wh ich the accuracy of 
any dynamical explanation should be measured l|Williams and Martined 1200(1) . 

Two major mechanisms for shaping food webs have been considered. One is the selection of 

those food webs that support stable population dynamics, because the degree of stability 

observed in nature is hard to achieve with simple, randomly assembled food-web models l|Mavl 
The other mechanism is evolutionary dynamics, which is suggested by the 
observation that t he similaritie s in th e trophic role of species and their phylogenetic kinship are 
highly correlated l|Cattin et alll2004[) . Even though there has been plenty of work on dynamical 
food-web models incor porating population dynami cs, evolutionary dynamics, or combinations of 
both (for a review, see lDrossel and McKaneLl2003f) . non of these models could so far be shown to 
reproduce fo od- web structure to the accura cy of the more recent descriptive models, namely the 
niche model ijWilliams and Martined 1200(1) and the nested hierarchy model ijCattin et all 12004^ . 
Purely evolutionary models (e.g.. lAmaral and Meveri 1999) have been developed mainly for 
understanding the dynamics of the total number of species, and not the detailed food-web 
structure. Work on m o dels combining evolu t ionary and pop ulation dynamics [e.g. 
ICaldarelli etHI lll998l) : iDrossel et all ll200ll) : lYoshidal jjoog)] has not yet arrived at conclusive 
results, one problem being high model complexity. A rec ent systematic tests of a plain 
population dynamical model bv iMontova and Solel (|2003) turned out clearly negative. 

Here we introduce a new model (called speciation model below) for the evolutionary dynamics 
of food webs, and show that its accuracy is comparable to the descriptive mod els. For an 
in-depth mathematical analysis of the model we refer to iRossberg et all 1120051) . 



2 Model definition 

We model the evolutionary dynamics of an abstract pool of species belonging to a given habitat, 
and the topology of the network of feeding relations between them. In order to separate the plain 
macroevolutionary dynamics from the complex conditions that determine fitness, evolution is 
modeled as being undirected. This simplification is here not meant to be a statement about 
nature, but is used as a technique to isolate the effect of particular processes. Population 
dynamics, for example, does not enter the model. 

The speciation model is most naturally described in terms of a continuous-time stochastic 
process. (See Appendix lAl for a description of our computational implementation.) Its dynamics 
consists of two parts, the dynamics of the community and the dynamics of the food web 
connecting the species. Table lists the model parameters. 
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2.1 Evolution of the species pool 

Only three kinds of macro-evolutionary processes are taken into account: speciations and 
extinctions within a habitat and adaptations of new species to the habitat. Anagenesis (evolution 
without speciation) is omitted for simplicity. Since changes in the environment can be fast on 
evolutionary time scales, it is for this model more appropriate to conceive a habitat as a 
particular set of environmental conditions supporting an ecological community, rather than just a 
specific location. 

Species at higher troph ic levels tend to be larger than at lower levels ijMemmott et all l2000l 
iLeaper and Huxhani [20021. Due to this correlation, a possible correlation between body size and 
speciation- and extinction rates leads to correlations be tween evolutio n rates and trop hic ley el. 
Arguments in favor of both, evolution rates increasing and decreasing l[Busbl . 

1 1993fl with body size have been put forward. Empirically, the question seems to be undecided. 
Generally one would assume that evolution rates are on the average proportional to some power 
M a of the adult body mass M. We note that a = is just as difficult to motivate as any other 
choice. Based o n the comparison of the observed distributions of the number of links to resources 
and consumers l|Camacho et all l2002bl). it can be argued that a < is more plausible in the 
context of the model proposed here l[Rossberg et all l200ljl . 

We associate each species i with a "speed parameter" Sj (triangles in Fig. ^| that determines 
the magnitude of its evolutionary rates ~ exp(si). The community dynamics is characterized by 
the values of the rate constants for new adaptations r\ , speciations r + , and extinctions r_ , the 
range [0, R] of allowed values for s, and the evolutionary speed/size dispersion D. The 
complexity of the processes driving evolution is accounted for by modeling extinctions, 
speciations, and adaptations of new species to the habitat as independent random events: For 
any short time interval [t, t + dt], the probability that a species i becomes extinct is r_ exp(si) dt 
and the probability that it speciates is r + exp(si) dt. When i speciates, a new species j with 
slightly different speed parameter sj = Sj + S is introduced into the community, where 6 is a 
zero-mean Gaussian random variable with vaiS = D. If + 5 exceeds the range [0, R], 
j = — (sj + S) or Sj = 2R — (si + 5) are used instead. The probability that a new species with a 



s 



speed parameter in a small range [s, s + ds] adapts to the habitat and joins the community is 
?*i exp(s) ds dt for all s between and R. 

2.2 Dynamics of the trophic network 

The food web is described by the topology of the network of directed trophic links connecting 
species. Parameters determining its dyna mics are the loopine s s A, the raw connectivity Cn, a nd 
the reconnection probab ility (3. Following lOohen et all ill 99(t . IWilhaTn7 and Martinezl (2000), 
and lGattin et all 1 2004ft . we take the fact into account that consumers are typically larger than 
their resource species. Empirical dat a show that they are, parasites excluded, at least n ot much 
smaller in terms of adult body mass l|Memmott et all l2000l ILeaper and Huxhaml l2002h . We 
assume that there is an upper limit h for the ratio of the adult body sizes of resource and 
consumer species. Resource species are at most h times larger than their consumers. Making use 
of the assumption that evolution rates ~ exp(s) scale as M a with a < 0, the difference of the 
speed parameter s of a resource and its consumer has an upper limit approximately given by 
A := \a\ x ln/i. The degree by which the trophic ordering is violated within the species pool is 
expressed by the loopiness parameter A = A/R. Summarizing, size constraints restrict the set of 
possible consumers of a species i to the species I with si < Sj + XR and the possible resources to 
the species h with Sh > Si — XR where < A < 1 (shaded area in Fig. GJ. Thus, cannibalism and 
trophic loops l|Polisl fl99ljl can occur. 

When a new species k is adapting to the habitat, a possible consumer of k becomes its 
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consumer with probability Co and, with the same probability, a possible resource of k becomes its 
resource. When a species i speciates (Fig.^l, the connections of the descendant species j are 
determined in three steps. First, all possible resources and consumers of i that are also possible 
resources and consumers of j are linked to j . Possible resources and consumers of j that are not 
possible consumers or resources of i are linked with probability Cq. Then, for each possible 
consumer or resource of j, the information whether it is connected is discarded with probability 
(3. Finally, a link is established with probability Co to any consumer or resource for which the 
information was discarded. Species that are both possible resources and possible consumers of j 
are treated like two species: a consumer and a resource. This procedure ensures that the average 
ratio of possible links to realized links remains at Co , in line with the premise that evolution is 
essentially undirected. 

The steady-state properties of the model are independent of the init ial conditions. 
Approximations for some basic properties can be obtained analytically l|Rossberg et all 12005^) . 
For example, the expected density of species has a constant value r-\ / (r _ — r+) along t he s axis. 
Thus, the observed broad distributions of logarithmic body masses Ce.g. lFenchel [19931 are here 
modeled by a simple boxcar function. The number of species S is on the average 
S = Rri/{r- — r+) and, with L denoting the number or links, the average directed food-web 
connectivity defined as C = L/S 2 is approximately C « Co (1 + 2 A — A 2 )/2. The fraction of 
species that entered the habitat by speciation is r + /r_, the rest (1 — r + /r_) entered as new 
adaptations. 



3 Model verification 

The speciation model was compared with seven of the best available empirical webs (see Table [5J , 
all of which exclude parasites, in line with our model. Collecting food-web data is not at all an 
easy task, and the data is sometimes criticised for being incomplete or inconsistent. But it was 
for this kind of data that the regularities mentioned in Section ^ have been found, which suggests 
that the data does carries substantial information characteristic for food webs. 

A simple, vis ual comparison of the connection matrix of the food web of Little Rock Lake 
(jMartinezl Il99lj) with a model sample (Fig. [2J shows that the model reproduces the 
characteristic mixture of randomly scattered blocks and isolated links found in empirical data. 
The empirical web has several omnivorous species of high trophic level, which leads to a high 
density of links near the left edge of the connection matrix. Such isolated, dense blocks are only 
occasionally found in simulations. They are not typical for empirical webs either. 

A systematic model verification was done based on twelve commonly used food-web statistics 
(see Appendix lB|l . We kept the parameters R = lnlO 4 and D — 0.0025 fixed, since their precise 
value has little effect (doubling either D or R changes the x 2 values computed below by less than 
one.), and set r_ = 1 without loss of generality. Depending on r + , A, and (3 (see below), 
parameters r\ and Co were chosen so that estimates of S and C after data standardization 
(Appendix El match with the species abundance S e and connectivity C e of the empirical webs. 
Maximum likelihood estimates (see Appendix Q assuming a multivariate Gaussian distribution 
for the twelve food- web statistics were used for the remaining model parameters r+, A, and (3 
(Table El- As a benchmark, fo od- web statistics were also computed for the niche mode l 
|williams and Martinez! Hp00) and for the nested hierarchy model ijCattin et all l2004|) . Neither 
has free parameters. As is shown in Fig. El there is good overall agreement between empirical and 
model data for all models. 

As a systematic measure for the goodness of fit, we computed the x 2 statistics given by 

X 2 = (Ve-v) T C- 1 (v c -v), (1) 
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where v e denotes the vector of the twelve food-web statistics for an empirical web, and v and C 
are the corresponding model mean and covariance matrix (Appendix Q . In terms of x 2 statistics 
(Table |SJ the speciation model (x§) is more accurate than the niche (xn) anc ^ the nes ted 
hierarchy model (xh)j even when differences in the number of statistical degrees of freedom 
(DOF = number components of v minus number of fitting parameters) are taken into account. 
For example, the overall p value for either niche or nested hierarchy model is less than 10 -24 
times that of the speciation model. In all but one case, the speciation model yields the lowest \ 2 - 
This is partially due to inherently larger model variances and partially due to a better accuracy 
of the model averages: in 46 out of 84 cases the average of the speciation model is closest to the 
empirical value (p < 0.001). Smaller webs are modeled particularly well. For the five food webs 
with S e up to 35 (after species lumping) all Xs values lie within the 5% confidence interval 
Xg < 16.9; and so does their sum X^Xs = ^0-5 < 61.7. 



4 Discussion 



In order to exclude the fallacies that the surprisingly good quantitative fit of the speciation 
model (1) does not actually depend on evolutionary processes or (2) is merely a consequence of 
the data standardization and evaluation procedures deployed, the model was also evaluated with 
r + = fixed, i.e. with all speciation processes inhibited. In this case the parameter (3 also does 
not have any effect, so that the only remaining fitting parameter is A. The model fits worsened 
considerably. We obtained \ 2 values ranging from 65 (Chesapeake Bay) to over 1000 (Ythian 
Estuary, Little Rock Lake). 

Of course this does not prove that there is no other explanation for the observations, but it 
shows that, for the good fit of the model, the speciation events are crucial. Two effects of the 
evolutionary dynamics are probably particularly important: (1) I t leads to approximate ly 
exponential distributions in the number of resources of a species llRossberg et all l2005|) . 
distributions that are similar to t hose obtained wi t h the nic he model by prescribing a beta 
distribution for the "niche width" ICamacho et alJ (|2002ala) . (2) Related species have similar sets 
of consumers and re sources, an effect that the algorithm of the nested hierarchy model 
l)Cattin et all l2004|) was set up to mimic. 

The speciation model is at least as accurate in its predictions as the niche model, w hich was 
itself a big improvement over the, by now historical, cascade model ijCohen et all fl99"o[) . The 
speciation model is the first th at involves a mechanistic explanation, confirming the earlier 
hypothesis IjCattin et all l2004f) that the structure of food webs can be unde rstood as the ef fect of 
a sequence of evolutionary events. Recent theoretical and empirical results llMcCannlEoorl show 
that the population dynamics in c omplex food w ebs can be stabilized by mechanisms such as 
weak links and adaptive foraging ( Kondohl l2003|) . This indicates that population-dynamical 
restrictions on food-web structure are less severe than one might first assume, making room for 
evolutionary dynamics to leave its traces. 
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Appendices 

A An Efficient Implementation of the Speciation Model 

An evaluation of the speciation model based on a straightforward discretization of speed 
parameters s and time t would require excessive amounts of computation time. Here the data 
structures and algorithms that were used in our, more efficient, implementation are described. 

A.l From Poisson processes to event sequences 

The main loop of our algorithm advances the dynamics of the food webs by one evolutionary 
event at each iteration. The question which event comes next is decided based on the following 
considerations: Call the quantity Oi := cxp(si) the raw evolution rate, or simply evolution rate of 
species i. We first note that the probability for new species with an evolution rate in the range 
[a, a + da] to adapt is 

ds 

ri exp(s) ds dt — r\o — — da dt = T\da dt. (A.l) 
da 



Thus, the total rate at which species adapt to the habitat is simply 

v\ ■= / rida = n (<r max - er min ) (A. 2) 



with CT m in := cxp(O) = 1 and <7 max := exp(i?). Since all possible events are assumed to be 
statistically independent, the probability that any new- adaptation, extinction, or speciation 
event occurs in the infinitesimal time interval [t, t + dt] is v tot dt with 



V t ot 



:= i/i + (r_ + r+) ^a u (A.3) 



where the sum runs over the indices of all species currently in the species pool. The time r to the 
next event follows the distribution 

P{t) = vtot exp(-f tot r). (A.4) 

If the probability for some event X to occur in an interval [t, t + dt] is vx dt, then the 
probability for event X to occur as the next event is 

— • (A.5) 

Vtot 

To see this, consider first the situation with only two kinds of events: event A occurring at rate 
va and event B at rate vb (both are Poisson processes). The joint probability distribution for 
the time ta to the next event A and the time tb to the next event B is 

Pab{ta,t b ) = va exp(-zv A TA) x v B exp(-^ B r B ). (A. 6) 

The probability that A occurs before B is 

/>oo /.oo 
P[TA <T B ]= dT A \ dT B PAB(TA,T B ) 

d,TAVAexp(-V A TA - v b t a ) (A. 7) 



-L 




VA 

v A + v B ' 
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The probability 1A.5|) is obtained by letting event A be event X and event B be any other event. 

For the sake of formal simplicity the range of allowed values for the speed parameter s is 
reduced from [0, R] to the half open interval [0, R[ without affecting numerical results. In order to 
decide which event comes next in the sequence of events, a random number p equally distributed 
in the range [0, vtot[ is generated. If p < v\, the next event is taken to be a new adaptation and in 
view of Eq. IjA.lJl the evolution rate of the new species is chosen as a random number equally 
distributed in the range [<7 m i n , cr max [- If v\ < p, the next event is taken to be either the extinction 
or a speciation of the species i that satisfies 

(r_ + r+) ^ < /> - < ^ + ( r - + r +) 2 (A - 8) 

Another random number is used to decide if the species i goes extinct [probability r_/(r_ + r+)] 
or speciates. 

A. 2 Data structures for keeping track of the network 

The connectivity between species is stored in a matrix of bits {bij} with bij = 1 is j eats i and 
bij = otherwise. In order to avoid having to re-arrange the full matrix bij with every event, the 
number of rows and columns of {bij} is chose much larger than the actual number of species is 
the pool. Indices are dynamically allocated with every new adaptation or speciation, and 
deallocated upon extinction. Unused indices are kept in a stack data-structure. The 
determination of the connectivity follows directly the model definition. 

A. 3 Sampling the steady state 

It is important to notice that the main loop of the algorithm, which simulates one evolutionary 
event per iteration, does not proceed with constant speed along the time axis t. Episodes with 
many species require more iterations per unit time than episodes with few species. In order to 
sample the steady state evenly in t, the an estimate t of t has to be used, since the true time 
between events is left undetermined with the method described above. Such an estimate can be 
obtained by starting with t — and incrementing t by the current expectation value of the time 
between events Ttot := ^tot a t eacn iteration. 

The portion of the food web that involves small species with large a reaches the steady state 
much faster than the portion that involves only large, slow species (small a). In the initial stage 
of the simulation, it is therefore sufficient to account for the portion involving only slow species. 
In practice, this is achieved by setting the effective value R' of the parameter R to a low value 
-Rstart (e.g., -Rstart = 1) at the beginning of the simulation and increasing R' after each iteration 
by a small amount proportional to r to t- When R' has reached R, the simulation is continued for 
some time, and then a steady-state food web is sampled. For the next steady-state sample, the 
complete food web is discarded, and the simulation is restarted with a R' = i? s tart- We started 
each run with an empty species pool. 

In order to check if the relaxation time is long enough to reach the steady state, we verified 
that, for the steady-state samples taken, the density of species along the s axis has the 
theoretically expected constant value n/(r_ — r+) for D = 0. 

B Food-web properties 

The food-webs properties that were used to characterize and compare empirical and mo del webs 
were: the clustering coefficient l)Camacho et all l2002bt iDorogovtsev and Mendesl l20f)2t) ( Clust in 
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Figure |3J); the fractions of c annibalistic species ^Williams and Martinez!. boOOf) (Cannib) and 
species without consumers l|Cohen et all Il99(i|) (T. top pred ators); the relative standard 
deviation i n the number of resource species jSchoenerf Il989(l (GenSD, generality s.d.) and 
consumers l|Schoeneil Il989|) ( VulSD , vuln erability s.d.); the web a verage of the maximum of a 
species' Jaccard similarity ijjaccardl Il908j) with any other species (jWilliams and Martinez! . 12000) 
[MxSim); the fraction of triples of species with two or more res ources, which have sets of 
resource s that cannot be or dered to be all contiguous on a l ine llCa^ in ei^ a ll 120041) (Ddiet); the 
average l)Cohen et all ll99o|) (aChnLg), standard deviation llMartinezl Il99ll) (aChnSD), and 
average per-species standard deviation l|Goldwasser and RoughgarderJ. Il993|) (a Omniv, 
omnivory) of the length of food chains, as well as the log 10 of their total number ijMartined . Il99lf) 
(aChnNo), with the prefix a indicating that these quantities were c omputed using the fast, 
"deterministic" Berger-Shor approximation (Ber ger and S hor. 1990J) of the maximum acyclic 
subgraph (MAS) of the food web. The number of non-cannibal trophic links not included in the 
MAS was measured as aLoop. When the output MAS of the Berger-Shor algorithm was not 
uniquely defined, the average over all possible outputs was used. 

The statistics were computed after standardization of food webs to allow a consistent 
comparison of data: Since in many records of empirical food webs the lowest trophic level is 
poorly resolved, all species without resource species were lumped to a single "trophic species" 
after dropping disco nnected species and before the usual lumping of trophically equivalent species 
llCohen et allll990|) . for both empirical and model data. 



C Maximum Likelihood Estimates, Goodness of Fit, and 
Confidence Intervals 

Maximum likelihood estimation of model parameters and a systematic characterization of the 
goodness of fit ideally require knowledge of the joint probability distribution of the food-web 
statistics in the steady state for any given set of model parameters. We used the methods 
described below to obtain satisfactory approximation in a computationally feasible way. 

Since species number S and directed connectivity C tend to fluctuate strongly in the steady 
state of the speciation model, it is impracticable to restrict sampling to those webs with S and C 
in close vicinity of the corresponding empiri c al val ues S c and C e as done by 
Wil liams and Martinez! l)2000(k ICattin et all l)2004j) . Instead, all webs with values of S and C 
differing by less than 30% from S e and C e were used for the statistical analysis. In order to 
estimate the steady-state mean v and covariance matrix C of the twelve food-web statistics 
conditional to S — S e , C — C c , we first computed steady-state mean u and covariance matrix D 
of the vector 

«=(;), (c.i, 

consisting of the twelve food web statistics v and the vector w = (S, C) T from N = 1000 samples 
within the ±30% range. Using w c = (S c , C C ) T and components of u and D _1 written in the 
block matrix forms 

«-(£)• °-' = (mZ mZ )• < C2 > 

estimates for v and C can be obtained as 

w e ) , (C.3a) 

(C.3b) 



v =u v + M vv M vw (u v 
C=M-L 
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The linear projection l|C.3|) is straightforwardly verified to be exact when u and, as a result, v 
have multivariate Gaussian distributions and N — > oo. Maximizing the corresponding likelihood 
function for the empirical statistics conditional to S — S e , C = C e , 



1 



V(2^) 12 |C| 



exp 



1 



(Ve-vfCrVc-v) 



is equivalent to minimizing 



X 2 + ln|C| 



(C.4) 



(C.5) 



with x 2 as given by Eq. in the main text. The maximum likelihood estimates for the model 
parameters r+, A, and were obta ined by minimizing l|C.5|) using the global optimization 
algorithm of ISchonlau et alJ l)l998j) . 

When the number of samples N used for estimating v and C in the formula for \ 2 (main text) 
is small, statistical errors of these estimates broaden the probability distribution of % 2 beyond 
the usual x 2 -distribution. In the simplest case % 2 (precisely iV^ 2 ) then follows Hotelling's T 2 
distribution, for which analytic expressions are known. We decided to use the infinite sample-size 
confidence intervals (in 95% of all cases \ 2 < 16.9 for 9 DOF and < 61.7 for 9x5 DOF) instead 
of the finite sample size intervals for N = 1000 (< 17.1 and < 65.1). This is the more conservative 
choice and, since our procedure for estimating v and C is different from the one assumed for 
Hotelling's T 2 distribution, it is not clear if these finite-sample-size corrections would still apply. 
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Table 1: Empirical food webs, fitting parameters, and goodness of fit 



Food Web 


key 


s c a 




r± 
r_ 


n 


A 


Co 


P 


AN 


AS 


AH 












number of 


statistical DOF: 


ll c 


9 


12 


Bridge Brook Lake 


BB 


15 


0.28 


0.91 


0.17 


0.12 


0.37 


0.059 


43 


13.9 


44 


(Havens. 1992) 
























Skipwith Pond 


Sk 


25 


0.32 


0.93 


0.21 


0.009 


0.53 


0.012 


83 


11.4 


64 


(Warren. 1989) 
























Coachella Desert 


Go 


27 


0.34 


0.96 


0.13 


0.006 


0.58 


0.014 


37 


3.9 


152 


(Polis. 1991) 
























Chesapeake Bay 


Ch 


27 


0.08 


0.96 


0.21 


0.25 


0.06 


0.029 


14 


7.5 


23 


(Baird and Ulanowicz, 1989) 




















St. Martin Island 


SM 


35 


0.14 


0.80 


0.92 


0.000 


0.23 


0.034 


17 


13.8 


12 


(Goldwasscr and Rouehgarden, 


1993) 


















Ythan Estuary 


Yth 


78 


0.06 


0.95 


0.67 


0.001 


0.08 


0.040 


58 


36.8 


82 


(Hall and Raffaclli. 1991) 






















Little Rock Lake 


LR 


80 


0.15 


0.99 


0.13 


0.025 


0.16 


0.006 


52 


24.2 


165 



(Martinez. 1991) 

a after species lumping (Appendix IB1 

b niche model ijWilliams and Martinezl [2000) 

c speciation model (this work) 

d nested hierarchy model ijCattin et all [2004) 

e the statistic Ddiet (Appendix El, which is always zero for the niche model, was excluded 



Table 2: Model parameters and their theoretical range 



model parameter 




range 


rate constant for new adaptations 
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n > o 


rate constant for speciations 


r+ 


r+ > 


rate constant for extinctions 


r_ 


r_ > r + 


total range of evolution rates 
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R>0 


rate/size dispersion constant 


D 


D > 


loopiness 
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< A < 1 


raw connectivity 
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< C < 1 


re-connection probability 
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[newSpeciation2 . eps] 



Figure 1: Illustration of the speciation model 

The food web is represented by a connection matrix (large square). Horizontal lines represent 
species as resources, vertical lines the same set of species as consumers. Trophic links are indicated 
by circles at the intersection points. The s-axis marks the evolution rate of each species. The 
process shown is a speciation. Species j evolves from species i. Most trophic links of j are copied 
from i, but some are modified. 
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Little Rock Lake 

x x 



Simulation Sample 



x . xxx . . x x 

x.xxx. .x x x x. . 

x.xxx. .x x x x. . 



. .XXXX XXX.X. . .X X. . 

. . XXXX X.X.X. . .X X. . 

. . XXXX X X.X.X. . .X. -X X. . 

t.XXXX X.X.X. . .X. -X X. . 



X.X.X. . .X. .X. . 



. .XXXX. XXXX 

. .XXX. .xxxxx 

. .XXXXXXXXX X. . 

. .XXXXXXXXX X. . 

. . . XX . . xxxxxx 

. .XXXXXXXXX X. . 

X. . . XXXX 

. .XXXXXXXXX X. . 

X. .xxxxx 

. .XXXXXXXXX x. . 

. .XX. XXXXXX X. . 



X . XXXX . XXXX . . 



.X. .X X. . .X. .X. . 



. .X. .X X.X.X. . 



. .X. .X X.X.X. . 



.X. .X.X.X X X. . 



. xxxx.x xxx. . 



. .X. .X X.X.X. . 



.X.X X.X.X. . 



,xxx 



XXXXXX. .X.X.X 

xxxxxx 

XX . X . . XX 



X X.X.X X.X X. . .XX. . .XX X X 

X...X..X X.X.X X X...X..X.XX X X 

.X..X X...X.. X.X.X X X X...XX.X.XX X X 

X. .X X.X XX. .XX X X X. . .X.X. . 

X X X.X XX. .XX X X X. . .X.X. . 

X X.X.X. .X.X X XX X X.X..X 

X. . .X X.X X.X. X.X X X 

.X..X X-.-X- .X.X.X X X X.X. X.X X X 

X...X..X...XX.X.X..X.X X X.X. X.X X X 

X...X..X X.X.X..X X.X X XX X X.X..X 



X.X. X.X X. . 



X. .X X. . .X. .X.X.X X X 

...X...X.X...XX.X.X.X.X.X...X X 

X. .XXXXX X X...X.. X.X.X X X 

...XX.X..XX.X XX...X.X...XX.X.X.X.X.X...X X XXX X X X. .X XI 

X.X. .XX X. . .X.X. . .XX. X.X. X.X.X. ..X X XXX X X X . . X XX 

. . .X. .XXXXX X X-.-X- .x.x.x x X X XXX X X X..X XX 

.XXX. . XXXXXX X X X XXX X X.X X. .X XX 

...XX.X..XX X...X.X...X..X.X.X.X.X...X....X XXX X X X..X XX 

. X . XX . X . . XX . X XX...X.X...X..X.X.X.X.X...X X XXX X X X..X XJ 

X.X. .XX X...X.X...X..X.X.X.X.X...X X XXX X X X. .X XX 

.XXXX.X. .XX.X XX...X.X...X..X.X.X.X.X...X X XX X..X X X.X. XX.. .XXXXX. XX. . 

. . .X. .XXXXX X X. . .X. .X.X.X X X X XX X..X X X.X. XX. . .XXXXX. XX. . 



.X..X X. . .X. .X.X.X X X X XX X..X X X.X.) 



. . . XXXXX . . 
. . XXXXX . . . 
. . .XXXXX. . 



.XX.X. XX X. . 



. .XX XX.X. . 

. . XX XX . X . . 

. . XX XX . X . . 



(.X.X. . -X. .X.X. . 
X.X. . -X. .X.X. . 



: . xxxxxx x..x..x.x..xx.. 



t.xx.x.xxxxxxxxxx. 



[samples . eps] 



Figure 2: Examples for connection matrices of food webs 

The connection matrix of the food webs of Little Rock Lake (|Martined . Il99l|) and a simulation 
sample, both after species lumping (see Appendix [B| for lumping procedure). Model parameters 
are as in Table [5] Little Rock. Each 'X' indicates that the species corresponding to the column 
eats the species corresponding to the row. The ordering of the species is such as to minimize the 
number of upper diagonal links and otherwise random. 
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Figure 3: Comparison of model and empirical data 

The graphs display the values of twelve food- web statistics (Appendix [B| of seven food webs (key 
in Table |5j) obtained from the speciatio n model (red stars), in compa rison with the empirical 



data (horizontal lines), the niche model l|yVilliams and Martinez . 



nested hierarchy model l|Cattin et al 
deviations. 



200(1) (blue squares), and the 



2004) (blue circles). Vertical lines indicate model standard 
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